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ABSTRACT 

We present the discovery of substellar-mass companions to three giant stars by the ongoing Penn 
State- Torun Planet Search (PTPS) conducted with the 9.2 m Hobby-Eberly Telescope. The most 
massive of the three stars, K2-giant HD 240237, has a 5.3 Mj minimum mass companion orbiting the 
star at a 746-day period. The KO-giant BD -f48 738 is orbited by a > 0.91 AIj planet which has a 
period of 393 days and shows a non-linear, long-term radial velocity trend that indicates a presence 
of another, more distant companion, which may have a substellar mass or be a low-mass star. The 
K2-giant HD 96127, has a > 4.0 Mj mass companion in a 647-day orbit around the star. The two 
K2-giants exhibit a significant RV noise that complicates the detection of low-amplitude, periodic 
variations in the data. If the noise component of the observed RV variations is due to solar-type 
oscillations, we show, using all the published data for the substellar companions to giants, that its 
amplitude is anti-correlated with stellar metallicity. 

Subject headings: planetary systems-stars: individual (HD 240237, BD +48 738, HD 96127); brown 
dwarfs: individual (BD +48 738) 



1. INTRODUCTION 

Searches for planets around giant stars offer an effec- 
tive way to extend studies of planetary system forma- 
tion and evo lution to st ellar masses substantially larger 
than 1 (iSato et al.l 2003. 2008a, 2010; Hatzcs ct ^ 
2006VD61hnger et al."2009; Nicdzicls ki et al.ll2007LI2009aF 
Liu et al.. 2009). Unlike their progenitors on the main- 
sequence, the cool atmospheres of evolved stars produce 
many narrow spectral lines, making a radial velocity 
(RV) measurement precision of < 10 m s~^ possible. 

There are now over 30 known planets and brown dwarf 
- mass companions detected around giants, providing 
statistics that are becoming sufficient to constrain the 
efficiency of planet formation as a function of mass and 
chemical composition of these evolved stars. For exam- 
ple, there is an ongoing debate as to whether the cor- 
relation between stellar metallic ity and planet frequency 
present in main-sequence stars (iFischer fc Valentill2005| ) 
is weaker or even entirely absent in giants, possibly due to 
the dilution of metal s as the convection zones deepen dur- 
ing giant evolution ("Pasqui ni et al]|2008t iZielihski et al.l 
[2011 in prep. I : iQuir rcnbach et al.ll2011[ ). A recent study 
of a homogeneous sample of subgiants by iJohnson et al.l 
(|2010) reveals that such correlation does exist in the case 
of these stars. 

Another interesting phenomenon is the observed de- 
ficiency of short-period planets around K-giants, with 
the closest one detecte d so far ha ving an orbital radius 
of 0.6 AU (Nicdziclski et al.| [2009a). This is in con- 
trast to subgiant systems, which include a few known 
"hot Jupiters" ((Johnson et al.ii201Q) . Models of the dy- 

^ Department of Astronomy and Astrophysics, the Pennsylva- 
nia State University, 525 Davey Laboratory, University Park, PA 
16802, sgettel@astro.psu.edu, alex@astro.psu.edu 

^ Center for Exoplanets and Habitable Worlds, the Pennsyl- 
vania State University, 525 Davey Laboratory, University Park, 
PA 16802 

^ Torun Center for Astronomy, Nicolaus Copernicus Uni- 
versity, ul. Gagarina 11, 87-100 Torun, Poland, An- 
drzej.Niedzielski@astri.uni.torun.pl 



namics of planetary systems around evolving giants sug- 
gest that the effect may be due to the tidal capture 
and engulfme nt of planets at the in itial orbital radii 
of 0.3-0.6 AU (ISato et al .ll2008bl: [yiU aver fc LiviQ,200l 
iNordhaus et al.ll2010[) . 

As a consequence of the observed correlation of 
the frequenc y of occurrence of massive planets with 
stellar mass (iLovis fc Mavori I2007t lOmiva et ahl 120111: 
IJohnson et al.ll2010D . there has been a growing number 
of det ections of brown dwa r f - mass companions to gi- 
ants (ILovis fc Mavoii 120071: iLiu et all l2008t ISato et al.l 
I2008at iDollinger et al.l I2OO90 . These include discover- 
ies of t wo-companion systems (jNiedzielski et all l2009bl : 
IQuirren bacli ct al. 2011), in which the orbiting bodies 
may be brown dwarfs, or super massive planets created 
via the core collapse process (|Ida fc LinI 120041 ). Re- 
cent studies have also uncovered several systems where 
a dwarf or subgiant is orbited both by a Jupiter-mass 
planet and a close b rown dwarf (e.g. iWinn "etaII[20T0l: 
IBenedict et al.|[2010D . In this paper, we report another 
intriguing case of a system consisting of a Jupiter-mass 
planet around the giant BD +48 738, and an outer com- 
panion, which may be another star, but the currently 
available constraints do not exclude a possibility that its 
mass falls in the brown dwarf range. 

The RV variations over a wide range of timescales and 
amplitudes observed in giant stars can be caused by a 
variety of phenomena, including orbital dynamics, pro- 
cesses that are intrinsic to the star, such as stellar os- 
cillations or spot activity coupled with stellar rotation, 
and other effects that may still await identification. For 
example, as shown by the CoRoT and Kepler photome- 
try of giants, both radial and no n-radial oscilla t ions are 
common m these stars (dc Rid der et al 1 120091: iHekkeri 
[20T0t IBaudin et alllMlt .Gilliland et al.ll201lD . In fact, 
p-mode oscillations have been identified in high-cadence 
RV m easurements of several of them (e.g. Barban et al] 
120041 [Hatzcs ct al. 2007; Hatzcs fc Zechn icister 200^. 
Studies bv lHekker et al.l (|2006t ) and lHekker et al.l (j2008f ) 
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demonstrate that one can make a statistical distinction 
between the stellar and non-stellar origin of RV varia- 
tions in giants by correlating them with stellar proper- 
ties, such as the B-V color and surface gravity. Here, 
we show two examples of such variability in the K2- 
giants, HD 240237 and HD 96127. Also, using the empiri- 
cally established dependencies betw een the effective tem - 
perature. Toff, and the B-V color (lAlonso et al}^199^. 
and T rff and surface gravity, log(g), (jHekker fc Mclcndcz 
[2001 . we investigate the intrinsic RV variability in a 
sample of K-giants with planets in terms of stellar metal- 
licity, [Fe/H], and find a suggestive anti-correlation be- 
tween [Fe/H and the observed stellar jitter. 

This paper is organized as follows. An outline of the 
observing procedure and a description of the basic prop- 
erties of the three stars are given in Section 2, followed 
by the analysis of RV measurements in Section 3. The 
accompanying analysis of rotation and stellar activity in- 
dicators is given in Section 4. Finally, our results are 
summarized and further discussed in Section 5. 

2. OBSERVATIONS AND STELLAR PARAMETERS 

Observation s were made w i th th e Hobby-Eberly Tele- 
scope (HET; IRamsev et al.l liggS) equipped with the 
High- Resolution S pectrograph(lTulL1998i') . in the queue- 
scheduled mode (iShetrone et al r i2007[ ). The spectro- 
graph was used in the R = 60,000 resolution mode 
with a gas cell (I2) inserted into the optical path, and 
it was fed with a 2" fiber. Details of our survey, 
the observing procedure, and d ata analysis have been 
described in detail elsewhere (jNiedzielski et al.l 120071 : 
iNiedzielski fc Wolszczad [20081 ). 

The atmospheric parameters of the three stars were 
determined as part of an extensive study of the 
FTPS targets described in the forthcoming paper by 
iZielinski et"al[(|2011 in prep. I . Zll). Here, we give a brief 
description of the methodology employed in that work. 

Measurements of Tes, logf g), and [Fe/H] were based o n 
the spectroscopic method of lTakeda et al.l (|2002l[2005aD . 
in which the Fel and FeH lines are analyzed under the 
LTE assumption. We used equivalent width measure- 
ments of 190 non-blended, strong Fel and 15 FeH lines 
per star which resulted in the mean uncertainties of 13 
K, 0.05, and 0.07 for T^g, log(g), and [Fe/H], respec- 
tively. Except for [Fe/H], these uncertainties represent 
numerical err ors intrinsic to the iterative procedures of 
iTakeda et al.l (2005a. b). The actual uncertainties are 
likely to be 2-3 times larger. Uncertainties of the [Fe/H] 
measurements were computed as standard deviations of 
the mean Fel and FeH abundances and are thus more 
realistic. 

Since the Hipparcos parallaxes of our stars are ei- 
ther unavailable or very uncertain, only rough estimates 
of their masses and radii could be made. Initial es- 
timates of stellar luminosities were derived from the 
My values resulting from the empirical calibrations of 
iStraizvs fc Kuriliend (|198lD . The intrinsic color indices, 
(B-V)o, and bolometric corrections, BCy , were calcu- 
lated from the empirical calibrations of lAlonso et al.l 
(|I999). 

Stellar masses and ages were estimated by comparing 
positions of the three stars in the [log(i/L0), \og{g), 
log(Tcff)l space with the t heoretical evolutionary tracks 
o fGirardi et al.. (,2000l) and lSalasnich et al. (,200M for the 



closest respective values of metallicity. Stellar mass esti- 
mates derived from this modeling process were then used 
to adopt the final values of luminosity. Stellar radii and 
the associated uncertainties were then computed using 
the spectroscopic values of Toff and adopted luminosities. 

Obviously, uncertainties affecting the position of a star 
in the parameter space determine the corresponding er- 
rors in mass estimates. As the values of Toff and log(g) 
are relatively well determined, the largest uncertainties 
are contributed by stellar luminosity or parallax. As a 
useful guidance, the mean uncertainty in stellar mass es- 
timates for over 300 Red Clump giants analyzed in Zll 
using the above approach amounts to about 0.3 Mq. 

The stellar rotation velociti es were estimated by means 
of the iBenz fc Mavoii (|1984f) cross-correlation method. 
The cross- correlation function s were computed as de- 
scribed by INowak et al.l ()2010[ ) with the template pro- 
files cleaned of the blended lines. Given the estimates 
of stellar radii, the derived rotation velocity limits were 
used to approximate the stellar rotation periods. The 
parameters of the three stars are summarized in Table [TJ 



3. MEASUREMENTS AND MODELING OF RADIAL 
VELOCITY VARIATIONS 

RVs were me asured using the standard I2 cell calibra- 
tion technique (jButler et al.| [T996) . A template spectrum 
was constructed from a high-resolution Fourier trans- 
form spectrometer (FTS) I2 spectrum and a high signal- 
to-noise stellar spectrum measured without the I2 cell. 
Doppler shifts were derived from least-squares fits of tem- 
plate spectra to stellar spectra with the imprinted I2 ab- 
sorption lines. Typical SNR for each epoch was ~200, 
as measured at the peak of the blaze function at 5936 A. 
The RV for each epoch was derived as a mean value of 391 
independent measurements from the 17 usable echelle or- 
ders, each divided into 23, 4-5 A blocks, with a typical, 
intrinsic uncertainty of 6-10 m s~^ at la level over all 
blocks (Nowak ct al. 2011, in preparation). This RV pre- 
cisio n levels made it quite sufficient to use the IStumpffl 
([198 Oi) algorithm to refer the measured RVs to the Solar 
System barycenter. 

The RV measurements of each star were modeled in 
terms of the standard, six-parameter Keplerian orbits, 
as shown in Figs. [1] [5| and U) Least-squares fits to the 
data wer e performed using the Levenberg-Marquardt al- 
gorithm ([Press et al.lll993 ). Errors in the best- fit orbital 
parameters were estimated from the parameter covari- 
ance matrix. These parameters are listed in Table [5] for 
each of the three stars. The estimated, ~7 m s~^ errors 
in RV measurements were evidently too small to account 
for the actually measured post-fit rms residuals (Table 
[2|). As outlined in Section 1, giant stars are subject to 
atmospheric fluctuations and radial and non-radial pul- 
sations that may manifest themselves as an excess RV 
variability. In particular, the p-mode oscillations, typi- 
cally occurring on the timescales of hours to days, are 
usually heavily under-sampled in the RV surveys of gi- 
ants and can account for a significant fraction of the ob- 
served post-fit RV noise ([Hekker et all [20061) . A rough 
estimate o f such variations can be mad e using the scaling 
relation of lKieldsen fc Beddind ([1995D , which relates the 
amplitude of p-mode oscillations to the mass, M, and 
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the kiminosity, L, of the star: 



-(23.4 ± 1.4) 



(1) 



where the aniphtude, Vqsc is in cm s~^. 

The statistical significance of each detection was as- 
sessed by calculating false alarm probabil ities (FAP ) us- 
ing the RV scrambling method (Wright et al .1120071 and 
references therein). FAPs for HD 240237 and HD 96127 
were calculated for the null hypothesis that the plane- 
tary signal can be adequately accounted for by noise. 
Two FAP values were calculated for BD -1-48 738, first 
for the null hypothesis that the long-term trend can be 
treated as linear, then for the null hypothesis that it can 
be modeled as a fraction of a long-period, circular orbit 
and noise. 

3.1. HD 240237 

RVs of HD 240237 (BD -1-57 2714, HIP 114840) are 
listed in Table [3] They were measured at 40 epochs over 
a period of 1930 days from July 2004 to October 2009. 
The SNR values ranged from 161 to 450. The exposure 
time was selected according to actual weather conditions 
and ranged between 184 and 900 s. The estimated mean 
RV uncertainty for this star was 7 m s""'^. 

Radial velocity variations of HD 240237 over a 5-year 
period are shown in Fig. [l] together with the best-fit 
model of a Keplerian orbit. The post- fit residuals are 
characterized by the rms value of 36 m s~^ and the FAP 
for this signal is < 0.01%. The observed RV variations 
point to an orbiting companion, which moves in a 746- 
day, eccentric orbit, with a semi-major axis of 1.9 AU, 
and has a minimum mass of TO2 sin i — 5.3 Mj for the 
assumed stellar mass of 1.7 M©. For randomly oriented 
orbital inclinations, there is a '^90% probability that the 
companion is planet-mass. 

The post-fit residuals for this model exhibit a large, 36 
m s^^ rms noise. Assuming that this RV jitter is due to 
solar-like oscillations, the Kjeldsen-Bedding relation pre- 
dicts an even larger amplitude of ^45 m s^^ for this star. 
To account for these variations, we have quadratically 
added 35 m s~^ to the RV measurement uncertainties 
before performing the least-squares fit of the orbit. 

3.2. BD +48 738 

RVs of BD +A8 738 (AG -1-49 313) were measured at 
54 epochs over a period of 2500 days from January 2004 
to November 2010 (Table|3]). The exposure times ranged 
between 420 and 1200 seconds depending on observing 
conditions. The SNR ratio ranged between 90 to 342. 
The estimated median RV uncertainty for this star was 
7 m s~^. 

The RV measurements of BD -1-48 738 over a 6-year 
period are shown in Fig. [21 They are characterized by a 
long-term upward trend, with the superimposed periodic 
variations. The best-fit Keplerian orbit model for the 
periodic component, with a circular orbit approximation 
for the observed trend, is given in Table [5] For the as- 
sumed stellar mass of 0.7 A^q, the planetary companion 
has a minimum mass of m2 sin i = 0.91 Mj, in a 393- 
day, moderately eccentric orbit with a semi-major axis 
of 0.95 AU. The FAP for this signal is < 0.01%. Within 



errors, this solution is not sensitive to a particular choice 
of eccentricity to model the outer orbit. 

A simultaneous characterization of the two orbits is 
difficult, as the outer companion's period is much longer 
than the time baseline of our measurements. However, 
the curvature of the observed trend is statistically signif- 
icant, with a FAP value of 0.01%. This makes it feasible 
to constrain a range of possible two-orbit models by fit- 
ting for parameters of the inner planet over a grid of 
fixed solutions for the long-period orbit. We have chosen 
to search for a minimum over a wide range of val- 
ue s in the a — m sin i — e space following the method 
of iWright et"all (|200l . The results are shown in Fig. [3] 
with contours marking the 1, 2 and 3-cr confidence levels 
computed under the assumption that the post-fit residual 
noise has a Gaussian probability distribution. 

To place limits on a range of the wide orbit models 
allowed by the available RV data, we considered solu- 
tions within the 2-a contour and further restricted them 
to those with e < 0.8, above which the occurrence of 
low mass stellar compa nions to F7-K primaries falls off 
(jHalbwachs et al.l 120051 ) . For the same reason, models 
involving plan et-mass companion s were limited to those 
with e < 0.6 (|Butler et al.ll2006f ). We tested the stabil- 
ity of these orbits, by carrying out many N-body sim - 
ulations using the MERCURY code ()Chambersl ll999D . 
Orbits from the grid in Fig. [3] were integrated for 10^ 
years over a range of a, msin i and e values within the 
2-a contour. 

As Fig. |3] clearly shows, the curvature of the long- 
period orbit is not yet sufficient to produce tight con- 
straints on the outer companion to BD -|-48 738. Al- 
though only a narrow range of massive planets, out to 
a ^10 AU, is still allowed, it is much more likely that 
the companion's minimum mass is in a brown dwarf or a 
low-mass star range and that its orbit is at least moder- 
ately eccentric. Continuing observations of the star will 
help placing tighter constraints on the orbiting bodies in 
this interesting system. 

3.3. HD 96127 

RVs of HD 96127 (BD +45 1892) were measured at 
50 epochs over a period of 1840 days from January 2004 
to February 2009 (Tabic E]). The SNR for these mea- 
surements ranged from 80 to 494. The exposure time 
was selected according to actual weather conditions and 
ranged between 72 and 1200 s. The estimated mean RV 
uncertainty for this star was 6 m s~^. 

The radial velocity variations of HD 96127 over a 5- 
year period are shown in Fig. |4l together with the best- 
fit model of a Keplerian orbit. As with HD 240237, the 
post- fit residuals show a large, 50 m s^^ rms RV jitter, 
comparable to the Kjeldsen-Bedding amplitude estimate 
of 45 m s~^. To account for these variations, we have 
quadratically added 45 m s""'^ to the RV measurement 
uncertainties. The best-fit parameters indicate a planet 
candidate with a minimum mass of 4 Af j, and a 647-day, 
1.42 AU, moderately eccentric orbit. The FAP for this 
signal is < 0.01%. 

4. STELLAR PHOTOMETRY AND LINE BISECTOR 
ANALYSIS 

Photometry and line bisector analysis, which are effi- 
cient stellar activity indicators, have been routinely used 
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to verify the authenti city of planet detec tions by means 
of the RV method (e.g. lQueloz et aLll200lD . In particular, 
when contemporaneous photometric data are not avail- 
able, it is important to correlate the time variability of 
line bisectors and RVs, in order to investigate a possible 
contribution of stellar activity to the observed RV be- 
havior. For example, even for slow rotators like K-giants 
discussed here, a stellar spot only a few percent in size 
coul d introduce li ne profile variations as large as ~ 100 m 
s~-^ (|Hatzesll2002| ). which would be easy to detect. 

Although a correlation between line bisector and RVs 
does provide evidence that the intrinsic stellar processes 
contribute part or all of the observed RV variability, 
many sources of jitter cannot be readily identified in a 
bisector analysis. This is because the bisector changes 
may be too small to measure distinctly from the appar- 
ent RV shift or they are indistinguishable from temporal 
changes in the spectrograph line spread function. Conse- 
quently, the lack of such a correlation is a necessary but 
not a sufficient condition for demonstrating that mea- 
sured RV variations are of a planetary origin. Additional 
information, s uch as a distinctly Kcplerian shape of the 
RV curve (Frink et all [2002: Zcchmcister et al. 2008) or 
a str ict persistence of the RV periodicity over many cy- 
cles ^Quirrenbach et aLll2011[ ). is necessary to fully verify 
such detections. 

In order to investigate a possible contribution of stellar 
jitter to the observed RV periodicities, we have examined 
the existing photometry data in search for any periodic 
light variations, and performed a thorough analysis of 
time variations in line bisector velocity span (BVS) and 
bisector curvature (BC) for each of the three stars dis- 
cussed in this paper. The BVSs and BCs were mea- 
sured using the cross-correlation method proposed by 
[Martinez Fiorc nzano (2005.) and applied to our data as 
described in iNowak et all (|2010t ). For each star and 
each spectrum used to measure RVs at all the observing 
epochs, cross-correlation functions were computed from 
~1000 line profiles with the I2 lines removed from the 
spectra. The time series for these parameters and the 
photometric data folded at the observed RV periods for 
the three stars are shown in Figs. [5l[6|and[71 

In addition, using the scatter seen in the photome- 
try data and the rotational velocities determined as de- 
scribed above, we have estimated the amplitude of RV 
and BVS variations due to a possible presence of spots 
on the stellar surface (Hatzcs 200?) . The calculated val- 
ues range from < 5 m s^^ to < 13 m s^^ and from < 17 
m~^ to < 27 m s~^, for the RV and the BVS amplitudes, 
respectively, which is much less than the variability ob- 
served in the three stars. 

4.1. HD 240237 

The longest time span photometry available for 
this star comes from 120 Hipparcos measurements 
(jPerrvman k ESAII1997I ). made between MJD 47884 and 
49042. These data give a mean magnitude of the star of 
V = 8.30 ± 0.02. There are also 15 epochs of photomet- 
ric observations of this star available from t he Northern 
Sky Variability Survey (IWozniak et al.ll2004l ). Neither of 
these data sets exhibits periodic brightness variations. 

The mean values of the BVS and the BC for this star 
are 136.6 ± 51.4 and 29.9 ± 47.3m s~^. No correlations 
between the RV variations and those of the BVS and BC 



were found with the respective correlation coefficients of 
r ^ -0.15 for the BVS and r = -0.06 for the BC. 

4.2. BD +48 738 

The only photometric database availab le for this star 
origin ates from the WASP measurements ([Pollacco et al.l 
I2OO6D . The WASP data are contemporaneous with our 
RV measurements and give a mean value of = 9.41 ± 
0.03, with no detectable periodic variations. 

We have also derived bisectors and line curvatures for 
all the RV measurement epochs, and obtained the mean 
values of BVS = 149.6±30.2 m s'^ and BC = 49.4±27.4 
m s~^. Time variations of these parameters are not cor- 
related with the observed RV variability and yield the 
corresponding correlation coefficients of r = 0.11 for the 
BVS and r — -0.08 for the BC time sequences, respec- 
tively. 

4.3. HD 96127 

Because this star is relatively bright, the only photom- 
etry available for it is fro m 100 Hipparcos measurements 
("Perr vman fc ESAI I1997). made between MJD 47892 and 
48933. The mean magnitude of the star calculated from 
these data is F = 7.43±0.01. In addition, the time series 
exhibits a marginal {2a), 25. 2±0. 05-day periodicity with 
the peak-to-peak amplitude of 0.02 mag. Analysis of 
the photometric time series divided into shorter sections 
shows that this oscillation persists over the entire ^3- 
year span of the data and, despite its low amplitude, it 
has a FAP < 0.001%. On the other hand, it is not present 
in our RV data for this star and it would be hard to envi- 
sion how could it be related to the observed 647-day RV 
periodicity. Moreover, the 4.3-day period corresponding 
to the frequency of maximum power of p-mode oscilla- 
tions predicted from lKieldsen fc Bedding! (|1995[ ) is obvi- 
ously much shorter than observed. Evidently, additional 
high-cadence photometric observations of HD 96127 are 
needed to clarify the nature of this 25-day periodicity. 

The mean values of the BVS and the BC for this star 
are 243.3±52.5 and 85.2±44.5 m s^"'^ and the correlations 
between the RV and the BVS and BC variations are not 
significant at r = 0.21 for the BVS and r = 0.17 for the 
BC, respectively. 

5. DISCUSSION 

In this paper, we report detections of periodic RV vari- 
ations in three K-giant stars from the list of targets mon- 
itored by the PTPS program. The accompanying analy- 
ses of the photometric data, as well as of the time vari- 
ability of line bisectors and bisector curvature indicate 
that the most likely origin of the observed periodicities 
is the Keplerian motion of substellar companions to these 
stars. 

The most interesting of the three stars is the case of 
the KO-giant, BD -^48 738, which, in addition to a > 0.91 
Mj planet moving in a ~1 AU, ~400-day orbit around 
the star, exhibits a long-term RV trend, which indicates 
the presence of another, more distant companion. The 
emerging curvature of the long-term RV trend is not yet 
sufficient to determine whether the object has a substel- 
lar mass, or is a low-mass star. 

In this context, it is interesting to note that there has 
been a growing number of detections of substellar com- 
panions to giants that have minimum masses in excess of 
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10 Mj, making them can didates for either brown dwarfs 
or supermassive planets (iLovis fc Mavoij 
1 l2008at iDolhnger et all 
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iLiu et ai] 



2009 ). In two 



[20081 ISato et al 
cases, double brown dwarf - mass companions have been 
discovered suggesting the possibility that they originated 
in the circumstellar disk, similar to giant exo planet s 
(|Niedzielski et al.ll2009bl:lQuirrenbach et al.ll20H"l ). If the 
outer companion to BD +48 738 proves to be substellar, 
we may have another interesting case of an inner Jupiter- 
mass planet and a more distant, brown dwarf - mass body 
orbiting the same star. 

In principle, such a system could form from a suffi- 
ciently massive protoplanetary di sk by means of th e stan- 
dard core accretion mechanism ()Ida fc LinI I2004D . with 
the outer companion having more time than the inner 
one to accumulate a brown dwarf like mass. A more 
exotic scenario could be envisioned, in which the inner 
planet forms in the standard manner, while the outer 
companion arises from a gravitational instability in the 
circumstellar disk at the time of the star formation (e.g. 
IKraus et al.l[20lll) . In any case, it is quite clear that this 
detection, together with the other ones mentioned above, 
further emphasizes the possibility that a clear distinction 
between giant planets and bro wn dwarfs may be diffi- 
cult t o make (see, for example, iQppenheimer fc HinklevI 
[20091) . 

The low signal-to-noise planet detections around the 
K2-giants, HD 240237 and HD 96127 reported in this pa- 
per further demonstrate that the intrinsic RV noise in gi- 
ants restricts the utility of the Doppler velocity method, 
especially in the case of long periods and large orbital 
radii, for which its sensitivity degrades f or dynamical 
reasons. This effect has been studied by iHekker et al.l 
(j2006( ). who have shown that the rms RV noise, any, of 
K-giants has a median value of ~20 m s"^ and it tends to 
increase toward later spectral types. Similarly, a negative 
correlation has b een found between anv a nd the stellar 
surfa ce gravity (jHatzes fc CochranI 119981 : IHekker et al.l 
[2008h . 

A short timescale, low amplitude RV noise commonly 
observed in giants appears to be t he manifestation of 
unde r sampled solar-like oscilla tions (jHekker et al.ll2006l 
l2008t iQuirrenbach et all 120 lit ). In particular, the re- 
cent photometric studies of K-giants by the CoRoT and 
Kepler mission s jBaudin et al. 2011; Ci ardi et al.ll2011l : 
iGilhland et"all 120111: IHekker et al.l I2O110 have shown a 
prominent presen ce of such oscillations in t hese stars. 
More specifically, iKieldsen fc Bedding (|1995f ) have pro- 
posed that the velocity amplitude associated with the 
p-mode oscillations scales as Vosc ^ L/M (Eqn. 1), or 
equivalently as T^^/g, where Toff, and g are the effec- 
tive temperature and the s urface gravity of the star. 
As shown, for example, by IHekker et all (|2006l l2008| l 
and iDoUinger et al.l (|2011| ) , observations are in a general 
agreement with these scalings. 

In order to compare the observed trend in the RV 
data for K-giants to the Kjeldscn & Bedding (1995) pre- 
diction more quantitatively, we have used the empirical 
dependencies of B-V color and log(t7) on r^ff for giants 
(lAlonso et al.l[l999l : [Hekker fc Melendezll2007n to gener- 
ate a semi-analytical function that relates a^y to the ob- 
served B-V for different stellar metallicities ([Fe/H]). The 
observed 30 values of grv have been taken from the pub- 



lished discoveries of substellar companions to K-giants in 
Table |6l as post-fit residuals after the removal of the Ke- 
plerian signal from the RV data. This sample s election 
has been motivated by the work of IHekker et al.l ([2008), 
which shows that amplitudes of RV periodicities in K- 
giants are generally independent of stellar parameters, 
whereas the RV noise left over after their removal tends 
to correlate with the B-V color and anti-correlate with 
log(g). Consequently, our approach practically ensures 
that the remaining short-term noise exceeding the instru- 
mental RV precision is intrinsic to the star, given the ob- 
served absence of planet s with orbital radii smaller than 
0.6 AU around K-giants ( Johnson et al.ll2007l : iSato et all 
I2008at iNiedzielski et al.l |2009aD. Moreover, as the RV 
precision and sampling patterns characterizing the on- 
going GK-giant surveys are quite similar, any selection 
effects that could influence the RV jitter measurements 
can be safely ignored. 

In Fig. m the measured values of urv as a function 
of B-V are compared with the Kjeldsen & Bedding re- 
lationship expressed in terms of these two observables 
as described above. The apparent increase of RV noise 
with the decreasing stellar metallicity becomes more ev- 
ident, if plotted out directly. The calculated Pearson 
correlation coefficient for this relationship is -0.71. The 
origin of this trend is most likely related to the fact that 
higher metallicity (opacity) of the star lowers its temper- 
ature, which decreases the amplitude of p-mode oscilla- 
tions, while lower metallicity has the opposite effect. 

One interesting consequence of the observed Ufiy- 
[Fe/H] dependence in K-giants is that it must affect 
the analyses of a correlation between the plane t fre- 
quen c y and metallicity of these stars (P asq uini et al.l 
mM. IHekker fc Melende"^ [20071 : IGhezzi et al.ii2010[ ). No 
clear agreement has been reached so far, as to whether 
the observe d correlation is weaker than in the case of 
dwarf stars ([Fischer fc Valentil 120051 ). or not. The effect 
unraveled here creates an apparent preference for higher 
metallicity stars to be more frequent hosts to orbiting 
planets, so that the published studies must be biased by 
this trend. 
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TABLE 1 
Stellar Parameters 



Parameter 


HD 240237 


BD +48 738 


HD 96127 


V 


8.19 


9.14 


7.43 


B-V 


1.682 ± 0.029 


1.246 ± 0.046 


1.503 ± 0.014 


Spectral Type 


K2 III 


KO III 


K2 III 




4361 ± 10 


4414 ± 15 


4152 ±23 


7r[mas] 






1.85 ± 0.89 


log(s) 


1.66 ± 0.05 


2.24 ± 0.06 


2.06 ± 0.09 


[Fe/H] 


-0.26 ± 0.07 


-0.20 ± 0.07 


-0.24 ± 0.10 


log(L*/LQ) 


2.49 ± 0.20 


1.69 ± 0.33 


2.86 ± 0.47 


M*/M0 


1.69 ± 0.42 


0.74 ± 0.39 


0.91 ± 0.25 


R*/R@ 


32 ± 1 


11 ± 1 


35 ± 17 


y™t[Km/s] 


<1 


<1 


<0.5 


Prot[d] 


>1010 


>505 


>1515 



Parameter 



TABLE 2 
Orbital Parameters 



HD 240237b BD -1-48 738b HD 96127b 



P [days] 
To [MJD] 
K [m/s] 
e 

[dcg] 
m2sin{i) [Mj] 
a [AU] 

Post-fit rms [m/s] 



745.7±13.8 
54292.0±28.3 
91.5±12.8 

0.4±0.1 
108.1±21.8 
5.3 
1.9 
1.1 
36.0 



392.6±5.5 
54457.2±28.5 
31.9±2.6 
0.2±0.1 
358.9±31.1 
0.91 
1.0 
1.1 
16.0 



647.3±16.8 
53969.4±31.0 
104.8±10.6 

0.3±0.1 
162.0±18.2 
4.0 
1.4 
1.3 
50.0 



HD 240237 



100 




-100 



-200 L 



100 



E 



-100 



53500 54000 54500 

Epoch (MJD) 



55000 



Fig. 1. — Top; Radial velocity measurements of HD 240237 (circles) and the best fit single planet Keplerian model (solid line). Bottom: 
The post-fit residuals for the single planet model. 
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TABLE 3 
Relative RVs of HD 240237 



Epoch [MJD] RV [m s"!] aRv [m 



ooloo. oy ( yo 


OO. ( 


'7 Q 
I .O 


ooo4d.4oUoz 


59.3 


/ .o 


00004. ZoiOU 


-o. / 


o.y 


Ov50Do.Uyzo4 


-Do.U 


Q 1 
O.l 


Ov5DD4.zoiyz 


-do.d 


o.l 


oooyo.44oiy 


1 OO A 


O.D 


oooy /.4oi /4 


161.8 


8.3 


ooyUU.44ilo 


123.0 


6.1 


OoyUo.4i4Uo 


113.6 


8.9 


00900.45211 


129.3 


7.0 


04Uzl.Z0oiy 


61.8 


4.5 


o4U41.UOo4o 


102.9 


6.2 


o4U0 / -Uzyoo 


1 Tcr o 
1 ( O.Z 


O.U 


o4Uo/.ioioi 


173.4 


8.1 


o4U (O.LZt oy 


1 DO O 


O.O 


o4uyz.iUozy 


147.2 


0.7 


04iiU.UoUoo 


194.7 


8.1 


o4z ( ; .4UoUz 


OO. 1 


/ .o 


O42o2.oyozy 


49.3 


5.3 


o4ool.iOio4 


-34.7 


6.6 


04oDU.lo4o0 


-29.2 


6.1 


O44oo.l4oz4 


54.9 


C A 

D.4 


"^441: U . 1 1 41: 1 


1 Si /I 

10.41: 


/t 7 
4. / 


044/4.UoUyD 






o40zi5.44oc5i 


of .XJ 


5.5 


O4ooz.4oy2y 


51.8 


4.7 


o4ooo.4zooo 


43.6 


1 o 


54639.38022 


119.8 


8.2 


54725.35386 


147.5 


10.3 


54743.31509 


94.9 


5.1 


54758.09993 


145.2 


7.8 


54777.22166 


129.9 


6.8 


54792.16502 


175.7 


7.2 


54812.12646 


106.2 


6.7 


54844.04772 


190.2 


8.4 


55049.45087 


56.2 


7.6 


55073.23259 


-23.0 


11.5 


55097.18015 


19.1 


7.2 


55118.24693 


7.6 


16.6 


55118.29165 


8.1 


7.1 
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TABLE 4 
Relative RVs of BD +48 738 



Epoch [MJD] RV [m s"!] ctrv [m s"!] 



53033. 65040 


-721.3 


11.6 


53337.60326 


-633.0 


8.0 


53604.86380 


-615.7 


7.4 


53680.90190 


-544.3 


6.6 


53681.65521 


-552.3 


7.8 


53970.86815 


-553.0 


5.9 


54006.77500 


-519.2 


6.3 


54022.73748 


-526.3 


6.7 


54034.69078 


-470.7 


5.6 


54048.65158 


-471.8 


6.8 


54071.79441 


-483.0 


3.8 


54086.56153 


-495.7 


5.0 


54101.71295 


-493.5 


4.5 


54121.67129 


-481.9 


8.9 


54121.68284 


-489.5 


6.6 


54134.63773 


-497.1 


9.2 


54149.60381 


-515.6 


8.8 


54159.58385 


-513.4 


5.9 


54162.57425 


-503.5 


6.6 


54341.85828 


-491.4 


7.6 


54365.78790 


-480.4 


6.3 


54381.96192 


-453.8 


5.7 


54399.68128 


-458.0 


6.6 


54423.63082 


-427.9 


5.6 


54440.79597 


-433.9 


5.4 


54475.71081 


-449.5 


4.8 


54507.61437 


-451.5 


5.9 


54687.89250 


-417.5 


7.9 


54713.82938 


-432.2 


7.4 


54729.77925 


-409.6 


8.3 


54744.74995 


-383.3 


11.7 


54759.92066 


-403.3 


6.5 


54774.68442 


-402.6 


8.0 


54786.85953 


-411.2 


7.6 


54794.61064 


-361.8 


6.8 


54812.59536 


-377.5 


7.1 


54821.55344 


-378.8 


7.0 


54857.66375 


-375.1 


9.3 


54883.60178 


-336.9 


6.0 


55048.92340 


-397.6 


4.9 


55073.84747 


-407.6 


6.0 


55102.77619 


-417.5 


3.4 


55172.57767 


-357.1 


8.3 


55198.71716 


-352.9 


4.9 


55222.66132 


-343.9 


7.1 


55246.58682 


-300.8 


7.5 


55246.61535 


-308.4 


8.5 


55249.61437 


-328.3 


6.2 


55444.83373 


-387.4 


4.5 


55468.74896 


-368.4 


7.1 


55493.71079 


-406.5 


5.5 


55493.71948 


-404.3 


6.3 


55498.68716 


-347.8 


5.1 


55513.63941 


-366.4 


6.4 
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TABLE 5 
Relative RVs of HD 96127 



Epoch [MJD] RV [m s"!] ctrv [m 
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-860.1 
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-992.4 
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TABLE 6 
Giant Stars with Known Planets 



Name 


B-V 


log(9) 


[Fo/H] 


Mi,/MQ 


R*/Rq 


rms [m s 


Rn 4-20 2A'S7 


1.25 


3.17 


-1.0 


2.8 


49 


60 


BD +48 7SS 


1.25 


2.24 


-0.2 


0.74 


11 


16 


HD 1 1 ()77 


0.93 


2.90 


-0.21 


1.91 


10 


29 


HD 1 SI SQ 


1.47 






4 


50 


42 


HD 1 7nQ2 


1.25 


3.0 


0.18 


2.3 


11 


16 




1.11 


2.10 


-0.15 


1.13 


10 


18 


HD 81 fiS8 


0.99 


2.22 


-0.36 


2.1 


13 


24 


HD 96127 


1.50 


2.06 


-0.24 


0.91 


20 


50 


HD 1(12272 


1.02 


3.07 


-0.26 


1.9 


10 


15 


HD 104985 


1.03 


2.62 


-0.35 


2.3 


11 


27 


HD 1 1 9445 


0.88 


2.40 


0.04 


3.9 


20 


13.7 


HD 1 S9S57 


1.20 


2.9 


-0.13 


1.35 


12 


14 


HD 145457 


1.04 


2.77 


-0.14 


1.9 


10 


10 


HD 1 7'^41 fi 


1.04 


2.48 


-0.22 


2 


13 


18.5 


HD isnsi4 


1.00 


2.98 


0.20 


2.6 


9 


13 


HF) 940910 


1.65 


1.55 


-0.18 


1.25 


11 


25 


HD 24r)2'^7 


1.68 


1.68 


-0.26 


1.69 


20 


34 


NOC 242S Nn S 


1.21 




0.14 


2.4 




18 


NnC 4S49 Nn 1 27 


1.46 




-0.12 


3.9 




13 


beta, Gem 


1.00 


2.68 


0.19 


1.7 


9 




gamma Leo A 


1.13 


1.59 


-0.49 


1.23 


40 


43 


epsilon Tau 


1.01 


2.64 


0.17 


2.7 


14 


10 


iota Dra 


1.18 


2.24 


0.03 


1.05 


13 


10 


xi Aql 


1.02 


2.66 


-0.21 


2.2 


12 


22 


4 UMa 


1.20 


1.8 


-0.25 


1.23 




30 


11 Com 


0.99 


2.5 


-0.35 


2.7 


19 


25.5 


11 UMi 


1.40 


1.60 


0.04 


1.8 


24 


28 


14 And 


1.03 


2.63 


-0.24 


2.2 


11 


20 


18 Del 


0.93 


2.82 


-0.05 


2.3 


9 


15 


42 Dra 


1.19 


1.71 


-0.46 


0.98 


22 


14 


81 Get 


1.02 


2.35 


-0.06 


2.4 


11 
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Fig. 2. — Top: Radial velocity measurements of BD +48 738 (circles) and the best fit model consisting of a circular approximation for 
the long-period orbit, and a 390-day planetary orbit (solid line). Bottom: The post-fit residuals for the aboye model. 
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Gettel et al. 
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Fig. 3. — Contours of and e in a - m sin i space of best-fit orbits to the RV data of BD +48 738 with in gray scale. The solid contours 
mark the levels at which increases by 1, 4 and 9 from the minimum, which corresponds to the la, 2a and 3a confidence intervals. The 
dashed contours approximately mark the levels of the eccentricity of 0.2, 0.4, 0.6 and 0.8, from the inner to the outer contour, respectively. 
Orbits interior to the grey contour are found to be stable for at least 10'' years. 



HD 96127 
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Fig. 4. — Top: Radial velocity measurements of HD 96127 (circles) and the best fit single planet Keplerian model (solid line). Bottom: 
The post-fit residuals for the single planet model. 
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HD 240237 
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Fig. 5. — HD 240237b data folded at the best-fit orbital period. Prom top to bottom: (a) Radial velocity measurements with the best fit 
model (b) Bisector Velocity Span (m/s) (c) Bisector Curvature (m/s) (d) Hipparcos photometry. 




Fig. 6. — BD -|-48 738b data folded at the best-fit orbital period. From top to bottom; (a) Radial velocity measurements with the best 
fit inner planet model and the outer orbit removed (b) Bisector Velocity Span (m/s) (c) Bisector Curvature (m/s) (d) WASP photometry. 
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HD 96127 
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Fig. 7. — HD 96127b data folded at the best-fit orbital period. Prom top to bottom: (a) Radial velocity measurements with the best fit 
model (b) Bisector Velocity Span (m/s) (c) Bisector Curvature (m/s) (d) Hipparcos photometry. 





Fig. 8. — Top: Post-fit rms RV error for Kcplcrian orbit models as a function of the B-V color for giant stars with published substellar 
companions (open circles), the stars discussed in this paper (filled circles), and stars with new, unpublished planet candidates from this 
survey (triangles). Theoretical predictions of the p-mode oscillation amplitude as a function of B-V, derived from a modified Kjeldsen- 
Bedding relationship, are shown as a set of curves for five fixed values of [Pe/H]. The shaded region indicates the range of RV precision of 
the different giant star planet surveys. Bottom: The same data plotted as a function of the observed [Pe/H]. Also shown are the theoretical 
curves for two fixed values of B-V. 



